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'\^ [ Abstract 

j^ ■ We discuss the longitudinal structure function in nuclear DIS at small x. We work 

within the framework of universal parton densities obtained in DGLAP analyses at 
NLO. We show that the nuclear effects on the longitudinal structure function closely 
follow those on the gluon distribution. The error analyses available from newest sets of 
nuclear PDFs also allow to propagate the uncertainties from present data. In this way, 
we evaluate the minimal sensitivity required in future experiments for this observable 
to improve the knowledge of the nuclear glue. We further discuss the uncertainties on 
the extraction of F2 off nuclear targets, introduced by the usual assumption that the 
ratio F1/F2 is independent of the nuclear size. We focus on the kinematical regions 
relevant for future lepton-ion colliders. 



1 Introduction 

Nuclear effects on the structure functions measured in deep inelastic scattering (DIS) ex- 
periments [Tl[2] offer valuable information for understanding the dynamics of partons in the 
nuclear environment. At small values of the Bjorken variable x, such measurements pro- 
vide a clean experimental setup for studying the behavior of QCD at high energies (see the 
review [3] and references therein). 

A usual framework for the determination of nuclear parton distribution functions (nPDFs) 
is that of global analysis using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evo- 
lution equations [1H6]. The procedure, pioneered in [714T0]. is similar to the ones performed 
for the proton case: obtaining the parameters for a set of parton distribution functions 
(PDFs) at an initial scale Qq which best reproduce some given sets of experimental data. 
The quality criterion is provided by a proper definition of a x^-function. As in the proton 
case, the goal of these analyses is twofold: on the one hand to check the degree of com- 
patibility of different data sets within an approach of universal PDFs evolved by DGLAP 
evolution equations; and, on the other hand, to provide a tool to compute cross sections for 
other processes in terms of a released set of PDFs for the different parton species. Given 
the fact that the factorization theorems in QCD [Ill[T2] are expected to be broken more 
easily for nuclei than for free protons, these checks are clearly of great importance in the 
phenomenological analyses of nuclear collider data. The latest sets of nuclear PDFs (nPDFs) 
obtained by these global fits are available at next-to-leading order (NLO) accuracy [T314T6] . 
with state-of-the art error analysis using the Hessian method [T^[T^. All these sets use data 
on nuclear DIS with charged leptons and Drell-Yan in proton-nucleus collisions. Data on 
inclusive pion production at high-pt has also been included in [15J providing extra constrains 
for the gluons without introducing tension among different data sets. Furthermore, this com- 
patibility of data within a universal set of nPDFs has also been checked in neutrino DIS on 
nuclei |17)H It is important to emphasize, however, that the role of the gluon distributions 
in all these checks is rather marginal (this can be seen in the corresponding error bars for 
gluons computed in p3l[T5] ) and further checks of this universality would be most welcome. 

The main caveat of these analyses is that an initial condition for Q^-evolution - not 
motivated from QCD but parametrized in a form flexible enough to reproduce the available 
experimental data - is required (for a different approach, see [TUH2T] and references therein). 
In this situation, the predictive power of the corresponding PDFs is reliable only in the 
region of x covered by experimental data: the extrapolations of both the central value and 
even the uncertainty bands outside this region remain linked to the functional form of the 
initial condition used in the analysis. For this reason, the results are not reliable for those 
parton flavors probed well outside the values of x constrained by the fitted data. This is 
particularly severe in the small-x domain and especially for gluons [21EIII5] which, on the 
other hand, dominate the cross sections at high energies. 

In DIS on nucleon targets, it is well known that an additional constrain on the gluon dis- 
tribution, on top of QCD scaling violations, comes from measuring the longitudinal structure 
function F^. This quantity has recently been extracted at HERA [22l|23] and the resulting 



^See, however, Ref. [TB] for contradictory results. 



impact on constraining the small-x evolution within the global DGLAP fits is currently un- 
der discussion [23] . On the other hand, due to the poor determination of the nuclear gluon 
distribution, the measurements of F^ on nuclear targets would be of great importance both 
for constraining the glue and for studying the nuclear dynamics at small x [25]. The existent 
experimental data on Fl are sparse and limited to a reduced kinematical region (see [26] 
and references therein). Actually, studies within perturbative QCD [27] and model calcu- 
lations [21] show that the corresponding nuclear effects closely follow those on the glue at 
small X. 

Furthermore, a knowledge of F^ is required in order to extract F2 from the measurements 
of the DIS cross section. In the nucleon case, the most recent combined HERA data [28] 
provide the cross section in the region y > 0.6 where the limited knowledge about Fl 
introduces a large uncertainty in the extraction of F2. In nuclear DIS, one usually assumes 
that the nuclear effects on both F2 and F^ are the same, i.e. that the ratio Fl/F2 is 
independent of the nuclear size and taken to be the same as that in the nucleon [29] . 

The purpose of this note is to analyze, within DGLAP approaches at NLO, the predicted 
nuclear effects on the longitudinal structure function and the uncertainties introduced in 
the extraction of F2 from the DIS cross section by our a priori lack of knowledge of the 
nuclear effects on Fl/F2. This will be done in Sections |2] and [3] respectively. We will focus 
on the kinematical regions relevant for future lepton-ion colliders [30|[3T]. We end with some 
conclusions. 

While our work aims at setting the uncertainties which would be acceptable by a global 
DGLAP fit, let us note that at small x and small-to-moderate Q^, fixed-order perturbative 
QGD which we employ here may not be suitable to describe the partonic structure of hadrons 
and nuclei, see [3] and references therein. Non-linear evolution equations which imply a 
saturation of partonic densities, offer another alternative. Nevertheless, the study of their 
impact on the longitudinal structure function has been done only for the proton [32], with 
the studies for nuclear targets performed within saturation-inspired models [331435] . We will 
make some qualitative comments on the expected results of saturation on the uncertainties 
introduced in the extraction of F2 from the DIS cross section at the end of Section [31 

2 Nuclear effects on Fl 

We define the nuclear ratio of function / (where / = F2,Fi,g, . . .) for nucleus of mass 
number A at momentum fraction x and squared virtuality Q^ as usually: 

where p stands for proton. At small x, experimental data indicate that Rp^ < 1, commonly 
referred to as nuclear shadowing. 

In what follows, we will work at NLO in the zero mass MS scheme, using CTEQ6.1M 
PDFs in the protoncl [36] and the corresponding nuclear ratios for the different nPDFs given 



^For HKN07, the MRST98 set of free proton PDFs are being used as their code gives absolute nPDFs 



in p^14T5| [T9]. With these sets of PDFs we then compute the corresponding values of the 
structure functions F2 and Fi (see the expressions in e.g. [3Z]) both in the nucleon and 
nuclear cases. 

In pQCD at NLO in the MS scheme, the longitudinal structure function Fl has the neat 
expression 



Fl{x,Q' 
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where fg denotes the gluon PDF, /^'s the corresponding quark PDFs, and e^ is the charge of 
a quark of flavor k. With gluon PDFs dominating at small x, the nuclear effects on F^ should 
follow those of the gluons. This is clearly seen in Fig. [H where Rp^ and R^ are plotted 
for Pb {A = 208). These results agree with the ones in [27] for those sets used in both 
analysis. Thus a measurement of Fl on nuclear targets offers the possibility of quantifying 
the nuclear effects on the gluon distribution at small x, which are essentially unconstrained 
in present-day analyses. 
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Figure 1: Results for Fi(Pb)/Fi and for the gluon density in EPS09 (left), and HKN07, 
nDS and FGSIO (right), for Q^ = 4 (top) and 100 (bottom) GeV^. The uncertainty bands 
and bars in EPS09 correspond to the propagation of the errors from the used data using a 
Hessian method in the global fit — see Ref. jTSj for details. The limits of the bands and 
bars for FGSIO correspond to two different model implementation in the calculation of the 
initial condition and the solid line is an average of these two values. 



"with this set as a baseline. 



3 Uncertainties in extracting F2 of the nucleus 

For neutral current DIS, the measurements are usually given in terms of the reduced cross 
section, which is a combination of structure functions. In the single photon exchange approx- 
imation and neglecting, at low and moderate Q^, the contribution from electroweak boson 
interaction, it reads 



2^NC 
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Y+F2 
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where a is the electromagnetic coupling constant and F+ = 1 + (1 — y)^. 

Thus the extraction of the structure function F2 from a^^ demands some knowledge 
on the ratio Fl/ F2. The uncertainties introduced by an inaccurate knowledge of this ratio 
become significant for large y (e.g. for y > 0.6, y'^/Y-^ > 0.31). With the experimental 
information about Fi in the proton currently being rather limited [22||23]. the most recent 
combined HERA data [28] provide the cross section, and not F2, in the region y > 0.6. 

In the nuclear case, the situation is more complicated by the additional nuclear effects on 
Fl/F2- The information about the nuclear F^ is sparse and limited at small 0.013 < x < 0.03 
to Q^ < 1.25 GeV^, moreover solely for rather light nuclei (see [26] and references therein). 
Previously, in the extraction of F2 from nuclear targets, the usual assumption has been that 
Fl/F2 is independent of the nuclear size and equal to that in the nucleoro. 

But it turns out that the recent DGLAP analyses [T3HTB| [TU] of nuclear parton densities 
give rise to nontrivial effects on the ratio Fl/F2 which, if neglected, may induce additional 
uncertainties on the extraction of the nuclear F2 from the reduced cross section. In order to 
estimate this uncertainty, we define the relative uncertainty 
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where F^ is the nuclear structure function extracted under the assumption of no nuclear 
effects on Fl/F2, while F^ is defined by Eq. (I3l). 

We show the results for Pb in Figs. |2]|11 We consider two kinematical situations, a 100 
GeV/nucleon proton or nucleus on a 20 GeV electron, and a 2750 GeV/nucleon proton or 
nucleus on a 50 GeV electron - which should roughly correspond to those collisions to be 
studied at the EIC [3^ and at the LHeC [3T] respectively. In Fig. ^ the uncertainty band 
corresponds to the one given by the application of the Hessian method to the nPDFs in 
the EPS09 parametrization [15], see that reference for details. In Fig. |3]and|l]we present 



■^E. g. using a parametrization 
and 0.6 < g2 < 20 GeV^. 



extracted from data on the proton and deuteron for 0.1 < x < 0.9 



the same quantity for the central values of the latest three NLO analyses of nuclear PDFs: 
EPS09 [TF|, nDS [13], HKN07 [13]. The trend of all these analyses is always similar with the 
different magnitude of the effect reflecting the different relative amount of gluon and quark 
shadowing for the three cases. Also shown in these figures are the results from the FGSIO 
parametrization [19]. As mentioned before, this set of nPDFs is built following a different 
procedure: the initial condition is obtained from the diffractive PDFs obtained in DIS with 
protons using the Gribov model of shadowing. Interestingly, this produces a correction with 
opposite sign to all other central values. 
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Figure 2: Results for the uncertainty in the extraction of F2 in EPS09 for LHeC (left) and 
EIC (right) kinematics. 



In order to understand qualitatively the numerical behaviour seen in Figs. [SHU we make 
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where we have assumed that F^ ^ F^ and the structure functions and nuclear ratios are to 
be evaluated at the virtuality corresponding to y = 1. Then, if the Q^-evolution of the gluon 
ratio is slower (faster) than that of the sea, the result turns out to be positive (negative) 
even at initial low scales and increases (decreases) initially with the scale but finally, at large 
Q^, it vanishes as shadowing dies out logarithmically with Q^. Examples of both behaviors 
can be seen in the figures. 
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Figure 3: Comparison of the F2 extraction-uncertainty in EPS09, nDS, HKN07 and FGSIO 
for the LHeC kinematics. For the latter set, the error bars reflect the variation resulting 
from the two options provided in there. 



This estimate also helps to understand the different behavior of FGSIO as compared 
with the other sets. According to Eq. ([6]), a gluon shadowing being much stronger than 
the corresponding one for quarks in the whole range of virtualities studied - as it is the 
case for FGSIO - would translate into a positive value of AF^. Notice that in a global fit, 
this quantity is also related with the logarithmic Q^-slope of Rp^{x,Q^) jlOlHT] which is 
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Figure 4: Comparison of the F2 extraction-uncertainty in EPS09, nDS, HKN07 and FGSIO 
for the ETC kinematics. For the latter set, the error bars reflect the variation resulting from 
the two options provided in there. 



constrained by existing data only for larger values of x than the ones studied here. 

On the other hand, the fact that for fixed x, A.F.f decreases with increasing Q^, is due to 
the kinematical factor y^/y+ in Eqs. (jlj) and ([5]), with y = Q'^/[xs), s being the Mandelstam 
variable of the lepton-hadron collision. Finally, the different behavior with energy of |AF^| 
for fixed Q^ and y, of the different sets (i.e. some of them show an increase, some of them 



a decrease), is due to specific features in tlie sets: tlie different slopes in x of the gluon and 
sea quark distributions, and tlieir beliavior at small x. 

While obtaining the analogous results in the framework of saturation models (see Section 
[1]) requires a dedicated study out of the scope of this work, some qualitative estimates can be 
done. Using the dipole model (see the corresponding expressions in e.g. [33]) and assuming 
scaling of the DIS cross sections with the variable Q'^ IQ\ with Q^ some squared saturation 
momentum which increases with increasing nuclear size and with decreasing x [121113] - 
as found in data and in non-linear QCD evolution -, the ratio of the longitudinal to the 
transverse DIS cross sections, ai/ar = (-Fl/-F2)/(1 — Fl/ F2), turns out to be generically 
a decreasing function of the scaling variable. Therefore, an increase in nuclear size (or a 
decrease in x) for fixed Q^ would imply a shift of the scaling variable towards smaller values 
and thus an increase in the ratiqj. 

All in all, the uncertainties introduced by the nuclear effects are sizable, rising up to ~ 10 
%, above all for small to moderate Q^ and small a; - as expected. This stresses the need of 
either measuring the longitudinal structure functions for nuclei or providing experimental 
results for the full DIS cross section in future experimental programs on lepton-nucleus 
colhsions [SUIEI]. 

4 Conclusions 

We have calculated nuclear ratios, defined in Eq. ([1]), with the most recent nPDFs for both 
nucleon and nuclei, shown the close correspondence between the nuclear effects on the glue 
and on Fj,, and found significant nuclear size dependence of the ratio F^ j F2 at small x. 
We went on to demonstrate how this theoretical uncertainty will effect the experimental 
extraction of the nuclear structure function F2- The resulting errors are largest, as high as 
~ 10%, in the most interesting kinematical region, namely at small x and moderate Q^ where 
future data can constrain the nuclear gluon distribution. This stresses the need for providing 
experimental data for nuclear DIS at future colliders in terms of reduced cross sections (total, 
charm and bottom), or, preferably, perform a collision energy scan to experimentally extract 
the nuclear longitudinal structure function. 
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